A series representation of the nonlinear equation for axisymmetrical fluid membrane 

shape 
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Whatever the fluid lipid vesicle is modeled as the spontaneous-curvature, bilayer-coupling, or the 
area-difference elasticity, and no matter whether a pulling axial force applied at the vesicle poles 
or not, a universal shape equation presents when the shape has both axisymmetry and up-down 
symmetry. This equation is a second order nonlinear ordinary differential equation about the sine 
sinip(r) of the angle ip(r) between the tangent of the contour and the radial axis r. However, 
analytically there is not a generally applicable method to solve it, while numerically the angle ip(0) 
can not be obtained unless by tricky extrapolation for r = is a singular point of the equation. 
We report an infinite series representation of the equation, in which the known solutions are some 
special cases, and a new family of shapes related to the membrane microtubule formation, in which 
simp(0) takes values from to n/2, is given. 

PACS number(s): 87.16.Dg,46.70.Hg,68.15.+e,02.40.Hw 
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Vesicles are closed surfaces of amphiphilic molecules 
dissolved in water, which form flexible bilayer membranes 
in order to minimize contact between the hydrocarbon 
chains of the lipid and water Recently, the vesi- 
cle shapes have attracted wide interest from different 
communities such as physics mathematics f|||U, 

chemistry and biology The naive model was 

given by Canham Q, in which only the surface bending 
elasticity is considered. An important progress was the 
introduction of the spontaneous curvature into the Can- 
ham's theory fiofl , and it was made by Helfrich with ana- 
logue to the spontaneous splay in the liquid crystal molec- 
ular layer. Since the membrane consists in two monolay- 
ers, if assuming that two layers depart from each other 
at a fixed distance 11 |l^], a so-called bilayer-coupling 
model was explored |13|. However, the individual mono- 



layers can in fact expand elastically under tensile stress, 
the so-called area- difference elasticity model Jl3|] was in- 
troduced and investigated p^ , |l5| ,p| . The latter three 
models are more realistic than the naive one, and have 
been intensively studied in recent years. A surprising 
common property of these three models is that they give 
exactly the same shapes [ p^l7|Ji~4] ] , while differing from 
each other in accounting for the shape transitions, such 
as budding and vesiculation transition, transition from 
prolates to oblates etc. Jlqjlq] . Furthermore, pulling or 
pushing the vesicle axially making it up-down and axis 
symmetrical, as it used in the experimental investigation 
of the tether formation Jl9|,[l7j , does not complicate the 
form of the mathematical equation for the general ax- 
isymmetrical vesicle in the spontaneous curvature model 
without force included. But the meanings of parameters 
need to re-specify. Except some numerical approach us- 
ing the powerful software Surface Evolver EfJ , usual nu- 



merical as well as analytical method is employed to solve 
the shape equation representing axisymmetrical vesicle. 
However, in analytical side [pT|,p2|| , there is not a general 
applicable method to study the equation. In numerical 
side, r = (cf. the following Eq.(l) ) is a singular point 
of the equation, so the shape around the point cannot 
be obtained unless by tricky extrapolation |16|j . In this 
paper, we are going to reveal a novel property of the 
equation: the equation has an equivalent series represen- 
tation, in which the two above problems do not present 
anymore. As a meaningful demonstration, a new family 
of shape related to the tether formation will be given. 

Since the detailed or key steps of the derivation of the 
same equation from different models is available in many 
papers, e.g. p7| , p3|j25| , we can start our further discus- 
sion from the equation itself with clearly specifying the 
physics meaning of each quantity in it. This may of 
course makes this paper more concise. And the equa- 
tion we start from reads: p^,0,p5[ 
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The quantities in this equation need some explanations. 
For an axisymmetrical shape, we can choose the symmet- 
rical axis to be the z-axis. The contour of the shape can 
be plotted in rz plane, with r being the radial coordinate 
denoting the distance from the symmetric z axis. Then 
the tangent angle tp(r) of the contour is measured clock- 
wise from r axis. Using s to denote the arclength along 
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the contour measured from the north pole of the shape 
and the <p the azimuthal angle, we have the following 
geometrical relations: 



|l = cosip{r) 
|f = -sinip(r) 
z(r) — z(0) = — J tanip(r)dr 

n = (simpcos<j), simpsinip, cosip) 



(2) 
(3) 



in which n denotes the normal of the surface. As to con- 
stants fc, C, Co, dp, and A, they have different meaning 
in different models. Taking the spontaneous curvature 
theory as example, fc the clastic modulus, C the inte- 
gral constant, cq the spontaneous curvature, Sp and A 
the two Lagrangian multipliers taking account for the 
constraints of constant volume and area, which may be 
physically understood as the osmotic pressure between 
the ambient and the internal environments, and the sur- 
face tension coefficient, respectively. When C = 0, Eq. 
(Ill) reduces to be the original Helfrich shape equation 
When c N, (2A + cg)/4 -> —L, Sp -► -6M, 
C — > — 4F, Eq.(f) becomes the equation representing the 
shapes with the same symmetry in area-difference elastic- 
ity model with the inclusion of the axially external force 
F (Eq. (19) in @). 

Although no one attempts to analytically solve it in 
a systematical way, we know that the following simple 
functions. 

sinip(r) = ar + b/r + c, and, ar + c^rlnr. (4) 

solve the equation p2| . Adjusting the constants a, b, and 
c in permissible ranges, the former form of the solution 
gives following shapes: 1) cylinder r = ro, sinip(ro) = 1; 
2) the axisymmetrical constant mean curvature surfaces 
( the so-called the Delaunay surfaces) with c = and an 
extension with c ^ 0; and 3) Clifford torus: simp(r) = 
—y/2 + r. The latter form of solution gives shapes from 
the self-intersecting oblates. oblate, sphere, prolate to 
very long capped prolate |26| . These analytical solutions 
not only have explained the existing but puzzling shapes 
|p|,p5|, but also had predicted something new. Clifford 
torus, for instance, was first theoretically predicted, then 
experimentally confirmed 

We will show that Eq.( 
for the physically interesting case, 
make a transformation for Eq. ([!]) 

ip(r) — arcsinf(r 



If) can be analytically solved 
For convenience, we 



and rescaling the unit such that k = 
Eq.(|I]) an equation for /(r): 

-2CV - 5pr 3 - 2/(r) - c 2 r 2 f{r) 
+2c rf 2 (r)+f(r) + 2krf'(r)~ 



(5) 

1 . Then we get from 

- 2Xr 2 f(r) 
2rf 2 (r)f'(r) 



+r 2 f{r)f 2 (r) + 2r 2 f"(r) - 2r 2 f 2 {r)f" (r) = (6) 

This is a second order nonlinear ordinary differential 
equation (ODE) not belonging to any well-studied math- 
ematical type. For our purpose to attack it, let us recall 



the Taylor series expansion of a function. As well-known, 
for any analytical function /(r) defining in a closed inter- 
val r e [ri, r 2 ], global Taylor expansion of f(r) around a 
point ro is 

, M ^ /W(r )(r-r ) fc , 



fc=0 



fc=0 



/( fc )(r )(r-r )* 
fc! 



(7) 



where R n (r - r ) = (r - r ) (n+1) / (n+1) (£) with n < 
£ < ra. Both the above series are exact and equivalent to 
each other. The first finite series is easy to serve the digit 
purpose, while the second infinite one is useful to obtain 
a closed expression once the general form is determined. 
To note that every axisymmetrical shape must intersect 
the equatorial plane at right angle tp — tt/2. This means 
that there exits a point r , such that /(ro) = 1. At this 
point ro, the two terms containing the second derivative 
in Eq. (||) cancel. We have then value of its first derivative 
/'(ro) of f(r) at point ro as 



/'(ro) = ± 



y/1 + 2Cr - 2c r + c 2 r 2 + 2Ar§ + SprJ 



(8) 



Since /(r) < /(r ) = 1, we have /'(ro) = nm (/W - 

l)/( r — r o) ^ if and only if ro is the maximum radius 
of vesicle. Since we treat f(r) a single valued function 
of r and ro is really the maximum radius of vesicle, only 
the positive /'(ro) is relevant. To note that the physi- 
cally interesting vesicles are either free or axially forced, 
and there is no force or torque acting on the waist, ft 
implies that the high rank derivatives of /(r) at ro as 
/' fe H r o))fc = 2,3... exist, and they can be obtained by 
the following method. Taking derivative with respect to 
variable r in both side, Eq.(^|) becomes a third order dif- 
ferential equation. Putting r = ro and substituting the 
value /(ro) (= 1) and /'(ro) (given by Eq. (g) ) into the 
third order equation, the third order terms cancel, and we 
can obtain the value of its second derivative /"(ro). Simi- 
larly, taking derivative with respect to variable r in both 
side of the third order derivative equation, we can ob- 
tain f^(ro), and so forth, all its higher rank derivatives 
/ (ro), k — 4, 5... can be determined. It is a remarkable 
result: the nonlinear ODE (||) can be transformed into 
an infinite series with the recurrence relation between co- 
efficients uniquely determined by the equation itself p8| . 
In the mathematical point of view, all its solutions with 
physical significance are obtained. 

In this paper, we only treat the simplest aspect of 
the solution without considering the congruence of the 
segments of shapes. Since the parameterization of the 
axisymmetrical surface is ip{r) rather than ip( s )i from 
Eq.(l) /(r) is a single value function of r in interval 
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4>{r) e [0,7r/2]. We can therefore expand /(r) in Tay- 
lor series directly. Conversely, once all Taylor coeffi- 
cients f^{ro)/k\,{k = 1,2,3,...), are known, the func- 
tion f(r) is determined from the relation (Q). In prin- 
ciple, we can give such a series for each shape. As 
checks of our method, we like to give two examples. 1) 
Cylinder with arbitrary radius r = ro. It is the case 
/(ro) = 1 provided all parameters k,Sp,X,co,C satisfy- 
ing 5prl + 2Xrl + k(c r - l) 2 = and C = @. It 
is the result discussed previously [^9 23 . 2) The De- 
launay surfaces f(r) = ar + b/r (a, b are two con- 
stants). These constant mean curvature surfaces ap- 
pear when A = Sp = C = 0. Now we have sim- 
ply / (?"o) = i( c o ~~ l/ r o) at position ro satisfying 
/(ro) = 1. We study the case with positive sign. We 
can easily find /W(r ) = (~l) (k+1) k\(c r - 2)/(2r§), 
(k = 2,3,4...). Then the solution can be exactly written 
as f(r) = ar + b/r with a = Co/2, b = (rp — cor 2 /2). It is 
the result discussed also previously ||22(| . One can verify 
this method using another analytical or numerical solu- 
tion as exercises, and can easily find that they are some 
special cases of the infinite series with special relation 
between the parameters. 

Now we wish to give new family of shapes. These 
shapes can be inferred from the works of three research 
groups in studying the meaning of the integral constant 
C. First by Zheng and Liu J23|, they noted that the con- 
stant C roughly meant that the constant pressure dif- 
ference Sp in Helfrich shape equation ]23| ] was replaced 
by Sp + 2kC/r 2 . Since that the all Helfrich shapes with 
axisymmetry are smooth at the poles r = |24|,[l6|], the 
constant C means a stress singularity at two poles. This 
singularity may lead to two horns at both poles. Sec- 
ondly, Jiilicher and Seifert claimed that the nonvanishing 
C connected to shapes with torus topology (30| . Thirdly, 
Podgornik, Svetina, and Zeks in |3l| ] and Bozic, Svetina, 
and Zeks in JljJ showed that a point axial force can lead 
to such a constant C 1 31 V\ . We have known that such a 
force can lead to microtubule formation of vesicle Jl7],[l9| , 
and in real experiment a finite force is enough to produce 
such structure. However, the expected shapes with two 
horns at both poles as exerting finite force have not been 
obtained theoretically. 

For our purpose, we treat a case in which the quantity 
in the square of Eq. (^) to form the square of a quantity 
with nonzero C. For our purpose, we choose the simplest 



parameters as C — cq, Sp — 0, A 



*/2, and the pos- 



itive sign of f'(ro). The first seven derivatives f^ k '{ro), 
(k = 0, 1,2, ..6) are 



Since the conformal invariance of the shape [lj,^,gG|.|17| . 
the magnitude of C does not matter. We can therefore 
choose C = 1. Then the Taylor series Eq.(^) can not 
give real and positive root of ro for equation /(0)= a 
negative number. It means that in this case the equation 
is not self-consistent. In contrast, it can give a unique 
and positive root for equation /(0)= a positive constant. 
To note that arcsinf(0) is the angle between the tangent 
of the contour at pole and the r axis. For a sequence of 
numbers of /(0) as 

N = {0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.972, 1}, 

(10) 

we have form the equation /(0) = N a corresponding 
sequence of unique and (semi-) positive roots for ro as 

r = {0,0.0095195,0.0189412,0.0282679,0.0375022, 
0.0466463, 0.0557024, 0.0646727, 0.0735589, 
0.0823627, 0.0886836, 0.0910858}. (11) 

Do not care the magnitude of these values. The energy 
of the vesicle depends on the shape due to the conformal 
invariance. For clearly specifying each shape, a conformal 
invariant, the so called relative volume v — V/ (4irR 3 /3) 
with R = ^/ A/Att is usually used (J], where V and A 
are the volume and area of the vesicle respectively. The 
one-to-one corresponding relative volumes are: 

v = {1,0.999971,0.999879,0.999725,0.999505, 
0.999213, 0.998844, 0.998386, 0.997819, 
0.997107,0.996449,0.996136} (12) 

In FIG. 1, the two limit shapes, sphere (thick solid 
line) and the prolate with the sharpest horn (thin line), 
are plotted. Other shapes with the horns of angles ip(0) 
satisfying < "0(0) < tt/2 will occupy the domain be- 
tween these two. 

For studying the relation between the "external force" 
C and the relative volume, we define a dimensionless 
external force as C(z — r ) with z being half of the 
length of the prolate, and rescale the relative volume as 
v' = 4.725(1 — v). We can find that both the dimen- 
sionless external force and the rescaled relative volume 
have nearly the same relationship to the zenith angle, as 
shown in FIG. 2. Results plotted in FIG.l and FIG.2 
clearly show that the zenith angle becomes sharper as 
the external force becomes larger, or equivalent, the less 
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the relative volume, the sharper the horn. 
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FIG. 1. The two limit shapes: sphere (thick solid line) and 
the prolate with right angle horns on two poles (thin solid 
line). They have relative volumes 1 and 0.996 respectively. 
Only a quarter of the contour is plotted. 



form simultaneously if no axially external force exerted. 
So, our study directly answers a question whether the 
constant C is related to the presence of the axially exter- 
nal force. This result is compatible with that obtained 
previously J3l],[l7],[30| . The obtained shapes are related to 
the formation of the membrane microtubule on the axi- 
ally strained vesicles, and finite force is enough to lead to 
such structure. Results also show that the shapes with 
less relative volume is easier to form such microtubule. 

Finally, we give two comments. First, even the Taylor 
series method is useful in obtaining and analyzing the 
solution to the equation, one should be careful to check 
the convergence of the series. In fact, we explore many 
cases with different parameter set of Sp, A, c , C, the con- 
vergence is sometime uncertain. Secondly, the presence 
of the constant C may not be the sufficient condition to 
relate the axial force Jl7],|2c|[il|, especially in the torus 
topology situation |3(J]. 
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FIG. 2. The correlation between the axial force C and the 
relative volume v. We plot the dimensionless force C(zo — ro) 
(solid line) with the zq being half of the length of the pro- 
late and the rescaled relative volume v' = 4.725(1 — v) (solid 
dots) . The larger axial external force and the smaller relative 
volume, the sharper angle of the corn. 

In summary, we have found that the second order non- 
linear different equation for axisymmetrical fluid vesicle 
shapes with up-down symmetry has an equivalent Tay- 
lor series representation, in which the recurrence relation 
of Taylor coefficients are uniquely and completely deter- 
mined be the equation. Using this representation, we can 
not only reproduce the known solutions, but also new 
ones. The new solutions discussed in this article give the 
shapes with horns at two poles. These horns can never 
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